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ABSTRACT 

A recent phenomenological model called the 
"Interacting Boson Model" (IBM) of collective levels 
in nuclei has been extremely successful in describing 
the rotational, vibrational as veil as transitional 
nuclei* The IBM treats the loulying collective 
states of a system of bosons possessing no intrinsic 
spin but able to occupy two levels with L=o and L=2. 

We have presented Hartree-Fock calculations, employing 
realistic interactions for several nuclei in the 
1 s-Od, Ip-Of and 1p-0f-0g shells with a view to examine 
some of the assumptions implicit in the formulation of 
the IBM. 
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CHAPTER I 


INTRODUCTION 


Nuclear Many -body Problem; Simple Shell Model and Configuration 
Mixing s- 

There is hardly a single worthwhile problem with 
realistic interactions in quantum mechanics which can be 
exactly solved for a many particle system. However approxi- 
mate solutions based on intuition have been of valuable 
help. An outstanding example of this is provided by the 
well-known atomic Hartree method. 

The success of the atomic Hartree method led the 

early nuclear physicists to believe that within a nucleus 

also an individual nucleon moves in some average potential 

( 11 ) . ^ 

produced by all other nucleons. In this way the shell 
model of the nucleus evolved and after that the shell model 
has gone through various stages of development* 

The average potential stated abo^ can be taken 
to be either an infinitely deep square-well or a harmonic 
oscillator potential. For both these potentials Schrodinger 1 s 
equation can be solved e xact ly and we gee various energy 
levels* In the single particle model , the nucleus is 
visualised to consist of filled shells consisting of 
maximum number of neutrons and protons permitted by the 
Pauli's exclusion principle as well as partially filled 
shells. This model is very simple in the sense that the 
nuclear properties are characterised by the last unpaired 
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nucleon. In an even-even nucleus all the nucleons are paired 

and so the spin of the nucleus is zero. In the case of odd- 

even nucleus the spin is given "by the last unpaired nucleon 

and in the case of odd-odd nucleus the spin is given by the 

unpaired neutron and proton. In the early stage of development 

of nuclear shell model one of the main objective was to 

explain the magic numbers. There are very strong experimental 

evidences that when a nucleus contains 2,8,20,50,82 or 126 

( 12 ) 

protons or neutrons a shell closure occurs. These numbers are 
called magic numbers . Let us see how the various stages of 
nuclear shell model could explain the magic numbers and other 
electromagnetic and nuclear properties of the ground state. 

The eigenstates are obtained in a spherically symmetric 
potential v(r) by solving the Schridinger ’ s equation, 


f v 


.m 


"e - v/ Y )} ^ ft/vO = +•- 


(T 0 


The solution of this equation is of the type 




/r, Y . 


In the case of infinite depth square well the radial 
wavefunction is given by the spherical Bessel functions which 


are regular at the origin 
So, in this case, 

U n i(r> = j 1 (kr) (I«3) 

From the boundary condition j-^CkR) = 0 where R is the 
radius of the square well, we get the energy eigenvalues. 

It is seen that the level 1 is 2(21+1) fold degenerate. 
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The closed shells occur for neutron and proton numbers 
2,2(21+1) with 1 = 0,1,2,.. i.e., at numbers, 2,8,18,20,3*+, 
k-0 ,58,68,70 , 92 , 106 , 1 12,1 38,156. So^ we see that t he infinite- 
depth square well can explain the magic numbers 2,8,20 but 
the magic numbers 28,50,82 and 126 are not reproduced by 
the square well* 

If we take harmonic oscillator potential then the 


radial wave function is given b; 


4 


12 J 




U n l<r) = 




, where. 


R* c w> 


N^xW- h> --W- 1 ) 




■(!■ 4) 


~ j-rx# C* 3 ) •" \ ^ which is the Lasuerre 

Polynomial and N is a normalisation constant. The energy 
levels are given by 

E h .~t3co ( bv, +- 1' - K) ( 1 , 5 ) 

n — - 1,2,3,** 1 = 0,1,2,.. 

The levels appear in the sequence of groups 1s;1p; Id, 2s; 
1f,2p; 1g,2d; 3s, etc. The closed shells occur at 
2, 8, 20, *+0,70, 112,168 

So, the harmonic oscillator potential also explains the 
magic numbers 2,8,20 but not the higher ones. 

Therefore we see that both the square well and the 
harmonic oscillator potential cannot explain the magic 
numbers beyond 20. To obtain the magic numbers 28,50,82, 
126, Mayer, Haxel Jensen and Suess suggested that a srpin 



orbit potential -f(r) (l.s) should be added to the centrally 

symmetric potential. This results in the splitting of 

3 = 1 ± i levels. ( s is the intrinsic spin of the nucleon) 

The choice of the sign in the potential is such that j=l + -j 

level is depressed relative to the j = 1 -g- level. This . 

explains all magic mumbers viz., 2,8,20,28,50,82,126. The 

infinite range harmonic oscillator potential can be brought 

closer to the finite square well if we cut off the range after 

a. finite value. This can be effectively be achieved by 

adding a ^ term to the potential. Here is negative. So, 

after this development the complete hamiltonian of a single 

necleon in a spherically symmetric potential is given by 

Tr 1- - - - ( 1G ) 

The spin-orbit splitting term has a relativistic orgin. 

actually a potential which is intermediate between the 

square-well and the harmonic oscillator well and which also 

has an experimental basis, is the Wood Saxon potential 

obtained by fitting the data on nucleon -nucleon scattering. 

This potential which is flat at the centre, and falls off 

(12J 

smoothly to zero at the edge of the nucleus is given by 
V(V ) ~ ~V r Q-t eXp/U- (-T-r)'I - ' - - . -(T '-?} 

where 0.5 x I0 _1 ^cm and R a-r 1 .33A sr x lo" 1 ^ cm, 

A being the mass number of the nucleus and 50-60 Mev. 

'.The single -particle model described above explains not 
only the magic numbers but the spin and parity of a large 



nurrb er of nuclei, provided we make the following assumptions. 

1 . Even-even nuclei, that is, nuclei with even Z and 
even II have total ground state angular momentum 1 = o + . 

2. An odd nucleus, that is, a nucleus with odd Z or odd 
II, will have a total ground state angular momentum and parity 
equal to the half-integral angular momentum j and the parity 
(-1)^ of the unpaired particle. 

3. An odd- odd nucleus will have a total angular momentum 
which is the vector sum of the odd neutron and odd proton 

J values. 

J = A + V 

The quantum number J is therefore an integer between the limit 



The parity will b e the product of the proton and neutron 

1+1 

parity that is, = (-1) n p . 

A partial justification for the above assumption arises 
from the tendency of the nucleon -nucleon interaction to give a 
J = 0 pair for two particles in (jm), (j-m) states. 

Although the simple, shell model predicts successfully 
the ground state spin and parity for many nuclei, it fails 
to give the ground state magnetic moment correctly. It 
can only give the limits between which the experimental depole 
moments lie. The two limits of the dipole moment when there 
is an unpaired nucleon in an odd-A nucleus, are 
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-rit ) 




f 


?6 + C Ti i)'h\ 


when J = H 
when J = 1--J 


where g Q and. g^ are the g-factors for the spin and orbital 
motion J " J 

The discrepancy between the experimental result and 
the calculated value of magnetic moment suggests that the 
simple shell model is not perfectly alright. The main 
drawback of the shell model is that there we have assumed 
a spherically symmetric potential of the nucleus, whereas, 
in reality most of the nuclei are deformed. 'The assumption 
of spherical symmetry is valid only in the immediate 
neighbourhood of a closed shell and when the qua.drupole 
moment is zero. In the case of a nucleus where the number 
of particles are such as to be in between two shells, the 
spherial symmetry is lost and the quadrupole moment has 
appreciable value. 

Due t o the incorrect assumption of spherical 

symmetry in the shell model the nuclear state is not given 

by the eigenstate of last unpaired nucleon, as the shell 
model* predicts. Rather, the eigenstate may be the linear 

combination of shell model eigenstates, and this eigenstate 

will represent a deformed nucleus. 

When there are a number of nucleons in a partially 
filled shell, outside a filled core, the nucleons can be 
in various configuration due to the residual interaction 
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between nucleon pairs. So the nuclear eigenstate will he 
the linear combination of the eigen st it es of various 
configurations for ?. particular total angular momentum. 

This is called conf iguration mixing . By considering 
configuration mixing we can calculate exactly the low-lying 
nuclear energy levels and the corresponding eigenfunctions. 

To illustrate the configuration mixing calculation, let us 
take the example of oxygen 18. It has two neutrons outside 
the closed core of 16. These two neutrons can b e in any of 

the three states Od^ 2 ? ^ s 'i/2’ 0d 3/2* ^ we v ‘ rant 

calculate the nuclear eigenstate corresponding to J = 0 we 
have to consider the case when the neutrons are in the 

i 2 ^ 

some j state. So the configurations will be \( od 5/2^/j = 0 

jci« , %2\ a nd !(0d^/«> 2 /-r ~ nuclear 

| US 1/2 ; /J = 0 an P U 3/2 ; / J = 0 

eigenstate will b e the linear combination of these states. 
The coefficients are obtained by diagonalising the matrix 
having matrix element of the type ) V I y where V is 

the two-body interaction. 

But, judging from the practical point of view , 
configuration mixing calculations can be done only for 
low-lying states of the nuclei and only when the configura- 
tion space is small. If we want to calculate the nuclear 
eigenstate for a high angular momenta or in a large confi- 
guration space ( i.e . , in the case of heavy nuclei) , the 
size of the matrix we have to diagonalise becomes enormously 
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large and prevents us to proceed further. We have however, 
a better alternative which is the deforme d Hart re e-Fock_ 
self-cons istent method. As yet the Hartree-Fock self- 
consistent met nod is known to he the most satisfactory 
approximate nuclear many -body technique for all sorts of 
nuclei where shell model calculations become impracticable. 
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CHAPTER II 


II.1 Hartree-Fock method as the best Zeroth-order treatment 
of the nuclear many -body problem s 


So, the Hartree-Fock method essentially is a many- 
body technique of finding out approximate nuclear eigenstates 
nowadays, two many body techniques are available in physics, 
broadly -speaking, (i) Canonical transformation approach 
(ii) Quantum field theoretic propagator approach. In 
canonical transformation approach the interacting system is 
transformed into a nearly nonintract ing system. In this 
way we can get two types of excitations. 

(i) Collective excitations (e.p., phonons in solid) 

(ii) Quasiparticle excitation (e.g., nuclei etc.) .In collective 

excitation we talk of the excitation of the system as a 

whole and in quasiparticle excitation we concentrate on a 

single particle surrounded by the cloud of other kind of 

particles. The quasiparticles ere almost non-interacting. 

We can naively say that the Hartree-Fock method is a kind 
of many -body technique which falls in the category of 

quasiparticle excitation • Here we talk about the 

Individual orbits but the orbits are noninteracting and 

deformed. The basic scheme is as follows? ' I 


'Y— /■ 

; - / ) 

’’ • 

Undeformed orbits 
with interaction 



- J 

( 

Deformed orbits with 
reduced, interaction 
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We go from the undeformed orbits with interaction 
to deformed orbits with reduced interaction, which is then 
neglected in the first approximation. Here, unlike the 
canonical transformation, we do not have to find out the 
eigenstate of deformed orbits by any transformation. We 
can immediately write down the hamiltonian of the deformed 
orbits since the deformation is quairupole type. The EF 
selfconsistent deformed orbits can also be interpreted., at 
a phenomenological level, in terms of the eigenstates 
of the IJils son hamiltonian 

jv— L C. -hM (./$ -f ' ~ • ' (]i. J-i ' 

The last term in the above expression contains the 
spherical harmonic yj which t-akes into account the 
quadrupole deformation. This hamiltonian is suitable 
when we deal with axial Hartree-Fock method, where the 
projection of the total angular momentum on the z-axis of 
the intrinsic frame of reference has a g ood quantum 
number. It should be mentioned here that deformations of 
most of the nuclear orbits are of quadrupole type. 

It is seen that when there is a number of nucleons 
outside a core of .filled, shells, the lowlying nuclear 
eigenstates are quasi -rotational t ype. The levels are 
not found to be exactly rotational in the sense that the’ 
energies of the various angular momentum levels are not 
exactly proportional t o J(J+1) but there Is an other term 
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which is proportional to J^(J+1)^ , but whose coefficient 
is small* actually, thi total nuclear eigenstate can be 
written as the product of intrinsic wavef onction and 
rotational wavefunction. This is called Bohr-Hottelson ' s 


unified model. In this model e.g., the state of spin <1 

belonging to the ground state rotational band of an even- 

( 1 ) 

even nucleus is described by the wavefunction 

i -™> njW- K: i" r ' ^ ‘(F 1 d 

The rotational motion is described by the symmetric top 


wavefunction 


fy (- 

i fs/f , 


1 1 ( f-rO ■ ■ 


IT. I ?) 


The intrinsic state i, <f:> t N / > is in the form of .antisymmetric 
si iter determinant containing the eigenstates of deformed 


orbits. The deformed orbits are given by Ae eigenstate 


of the Nilsson hamilt onian h 


Irw.hb = ebv> 

i \ / 


We’ll be dealing with the nucleons in the single 


major shell outside the closed core. This type of restrictive 
outlook has the justification that the energy difference 
between two major sheels is high enough to neglect the 
interaction between two major shells in comparison to the 
interaction between the nuclei within the same major shell. 

It is the task of Hartree-Fock method to determine the 
deformed orbits ) • This is done ^ - variational method 
given below. In this method we minimise the ground state 
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energy of the nucleus, which is the expectation value of 
the nuclear hamiltonian H in the Hartree-Fock states 


The nuclear hamiltonian H an:! the Hartree-Fock energy 

(m) — (jT 1 '■;) 


are given by 

i , _ ^ L *" 

n - - 7 ^ 

4.- i 


L L. 

a 


and 




A 


-c. 

> 


^ , f ' 1 


V 


T cah i>>+) "> 4-p- MV- > ■ is, 1 6 > 

c ' ' Z i 


respectively. 

if j V7 is the kinetic energy of the orbit y\ and 
|> / M ) =- ~ an a ntisymmetrised matrix element 

of the two-body interaction v(ij). 

It is advantageous to expand the orbits on a basis 


of known wave function. 

? v _ v { * \ 

~ / I t 

L~. '• I 


i-iv.) 


j 


The basis j may be e.g., the set of j t r. /■ of shell 
model states in a harmonic oscillator well. Any basis may be 
used to expand the orbits provided one is able to calculate 
the matrix elements of the interaction v with the basis 
states. The set of orbits j V? (both filled and unfilled) 
is assumed to form a set of orthonormal set of wave functions. 

i 


1 1 f *c 7 


b 


* 


) 


A ! ' 


Tv a' 


o - ■ 


(Il.i-s) 
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The wave function \ { 'fy is determined by the coefficient sC^ 
which become variational parameters. The energy will be 
stationary with the above normalisation condition when 



Here ^ is introduced as a Lagrange multiplier. With the 
help of the expansion of /\ „ the above equation may be 
expressed, in terms of known matrix elements of t and v in the 
known basis j. The derivs.tive may then be calculated 



This equation has the form of that of an eigenvalue po 
problem 



I 

where h is the hamilt onian given by the matrix elements 

A 


<i|hi7 -- <-, K 1-7 


/ 


r - < ' 


L 

A~ 


Hth' 


/ -f- 


/ 

** 




\ 

/ 


: a. h 


JJL. 1. !!.) 


h is the Hartree-Fock hamilt onian. The coefficients 
A Va are obtained f ran the above equation by a self- 

L r • tv 

consistent iterative method which is as follows j 
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(i) An initial set of values of the coefficients 
are taken . 

(ii) By these set of values the matrix elements of h are 

calculated from the en nation f T\ }. ill 

v”' ' 

(iii) The matrix is diagonalised and a new set of 

coefficients 1' /> are obtained. 

- j 

(iv) These ne' iT set of coefficients are utilised to 
calculate the matrix element and we hack to step (ii) 

In this way the iteration is continued until the 

coefficients in two successive iterations differ very 

little from e ich other. At every stage of iteration we 

can calculate also tho Hartree-Fock energy and compare 

its value in successive iterations. 

The choice of initial guess in Hartree-Fock iteration 

is very crucial because different types of initial guess 

lead to different Hartree-Fock minima. 

To get the rotational intrinsic Hartree-Fock 

state it is best to generate the initial guesses 

from the quadrupole operator i.e., we can generate the ^ 

.( 1 ) 1 

by ding on a.lis ins the matrix of quadrupole operator. 

From the axial Hartree-Fock method we can know the 
eigenstate of J z and the corresponding energy. That is, 
we can know some sort of average energy, and average 
wavef unction. Each value of total angular momentum projection 
corresponds to a rotational band. 



1.2 ungi 


il ar i I one n t urn Projection 1 


[ 1 , 7.18 ] 


Prora tho Kurtree-Foek ”-vef unct aon , 


; lie nuclei'.. 


ei S rr functions of GT,J Z ) can be obtained by J 
m one n t ura pro 3 ect ion__t e c nn ique_ . 


n~ular 


/, r\-^L V - - ' '(H 2 - 0 

l l — / ! ‘ 1/ U 


Let us 1 rite ^ — y 'T ' .* K 


where ^ . (v\P) is the Hartrec-Fock wave function for the 

K =0 bandV^Kere, we are expanding the average ^function 
obtained in Eartree-Fock method m terms of exgenst ate* 

J 2 and 

To find ; A cC s let us consider, 
i ( 

-nil-; 


=■ > Aj t ct -Ki-.'(f 
3 * *" 


{ 3 _ 2-2.; 


where <T ^ is a rotation operator around tho 2-axis 


and is tho d-matrix. p 

Thus, we get ^ 


is the angle of rotation. 


lus, . „-r 1 , x 


iA:iAR^L'(?)< Tk 

* -»r * * * 

7 > J 


TV 


7 - ( 1 - 2 . 3 ) 


iV 


Fr om this we can get A»s as 


x^u 


^ SP* o 


e 
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a technique to calculate <("'4',,, . .... 
is e.s follows. We cun construct 



K-o £ 




■ 1 tu X 


i A- / 


■ * V. J 


€ 


Fl 1, 


t K ../ 



/V 

1 *p\ 


unrotated 

determinant 


rotated 

determinant 


F )T example, let ^s consider the case of 2x2 determinant 
<? ’ 7 ■ € 


\ . (j) ^ o) 


Wt 




24 - . v C 4 


K ; 


4-< 0 

! ?„,y0) v4' 


■ (li- 2*6^ 


ch CO 

'iX = >/, 7 




■f h W 


A 




e 1 




--7 c>V t e ^. 

* H ^ 'U y K i * r » 
So, u - 


; a/ 




where, 

•i\4> !<H d * K/l 


<•1 7d’ /(rt<"k-> 


2 4<f 

4 K‘ 


4 


JU 







* l 


The energies corresponding to different J ? s cart be 


obtained as follows 


I" 




^TK 1 H I 



In equation - (11,2*^) we have obtained, 

i A d i F ^V-o i e'I Jv I W 

In a similar manner, we can set. 


{ fijd -Q-, 


£j\h) 


v E- 


2— „ a 


r^i. 


kro 


" fir i-9) 

get, h 

'^^hEk- o/^Kp)^- 


jVur t-d^HK rf ,)cE(p) d -n_ 


By this formula, the energy corresponding to each J in 
a band is determined. 


II. 3 The interacting boson model and investigation of 

the validity of its basic assumption in the frame- 
work of EF met hod: - 


it'uempts have been made at various stages to develope 

a truncation scheme for the large • c onfigur at ion space 

of shell model, so that shell model calculations become 

(8) 

practicable for medium and he,.vy-ma.ss nuclei. But such 
a truncation schemes have not yet been completely successful. 
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Recently a much simplified model has cone into picture, 
which can successfully reproduce the collect ive behaviour 
of nuclei, specially the medium .and heavy-mass nuclei. 

This model is called Interacting Bo s on iiode l. The central 
ir?ea of this model is this; The fermion system of the 
nucleons is grouper into pairs which rie bosons. The 
collective behaviour of the nuclei depends only on the 

' do, 

bosons having spins 0 and 2 i.e., the s-bosons and d-bosons* 

The s-bosons an:’ d-bosons interact among each other to form 

the multiparticle system. The Hamiltonians of interacting 

boson model are found in literature and such Hamiltonians 

have analytic solutions for energy spectra etc. But 

derivation of interacting boson model fh om shell model 

i.e., the microscopic description of the model has not 

yet been developed with much success. So, why the collective 

behaviour depends selectively on two particle states of 

110 ) 

angular momenta 0 and 2 only, Is not yet known. 

Hartree-Fock method was developed to do away with 
the difficulties connected with shell-model calculations, 
and here we do not need any truncation scheme for simpli — 
fication. But the basic assumption of interacting boson 
model can easily be judged In the frame-work of HF calcu- 
lations by switching off the two body interactions when 
their total angular momentum is other than 0 or 2 . We 
have found that 'the qualitative features of low-lying 
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spectra remain unchanged ’/hen we 50 from full interaction to 
selective interaction , involving two-particle states having 
angular momenta 0 and 2 only. Thus the basic assumption of 
the interacting boson model is justified. 
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CZ.J.PTH2 III 


C.iLCJL .TIOIUL ?d XlJURFS 


he hav" described th" Kurtree-Foch eruptions and 

'■.ng '.la.;- momentum pro 3 ec : ion in CP apt or II. Nov * r Q'll 

briefly describe the c.ilculational procedure ulongwith 

the L' echnicuee involved in computer programming. 

III.1 Storage of two-body matrix element s:- 

We have used matrix elements of Kuo -Brown realistic 

( 22 ) 

interaction available in liter -fcure . The two-body matrix 


elements of the form <^3^3, 


V ji : ' re stored in a 

three-dimensional array in the progr an, for easy reference. 
There is a special indexing procedure for this. Let us 
illustrate this indexing procedure taking the example of 
the s-d shell. In the s-d shell space, (0d^y 2 , 1 s -i/2 , 0<3 3/2^ 
We have three angular momenta which we can index as 
1 ,2,?. how, 3 1 ,3 2 on the bra-side of the matrix element 
can be indexed as 

a, = 3(^-1) + j 2 £111.1 .1) 

when 3^ and ^ vsxles from 1 to 3, 
lXj varies from 1 to 9* 

Similarly on the ^ e t s i de can de i nd e xed a s 

n 2 = 3C3-I )+ (IIX.1.2) 

For JT, we have n^ = 2J +1 (III.1 .3) 
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when J = 0 , T = 1 , n^ = 1 

J = 1 I=o n^=2 and so on. 

In s-d shell n^ can run from 1 to 9, from 1 to 9 and 
from 1 to lo. 

By this indexing procedure , also it is easy to incorporate 
the symmetry properties of the matrix elements. The 
matrix-elements liave the f ollowing symmetry properties. 

| Y 1 hy <JT = \^1 ^2 S Y ! ^3^V' JT 

j +j +J+T 
x (-1) 1 2 

^3^ \ 7 j 3 2^1/ IT = ^2^1 | Y I 3 3V>JT 


<^M2 | Y \ V 3 / jt = <^i 3 2| y |^y>jr 





I x ( 

_ i) 3 3 + V J+T 

9 V3 . 

V 

H \ _ 
^1^2 y jt 


7 j JT 

>< j | 

y2 3 i | 

V 

! M 3 / JT = 

<V 2 

7 1 ^\y jt 



X (-1) 

3 1 +;)2 + 33 + V 

\ Uhl 

V 

X 

, hhy' jt = 

\ 3 2 3 1 \ 

j Y J \^3y jt 


I 

i 

i 

i 


(III.1 .3) 


Xhi* u;u. 


nPUR 

iBRAfflf 


Bte. 
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III. 2 Definition of the basis 

The oasis In which the deformed orbital of the nucleon 
is expanded should be such thvfc the expansion converges and 
the matrix elements of two-body interaction in such basis 
are readily available. The harmonic oscillator wavefunction 
shows good convergence and in shell model basis the 
matrix elements are available. We have chosen the basis such 
that they are suitable for non-axial as well as axial 
Hartree-Fock calculation. 

In axial Hartree-Fock calculation, the cfe formed orbilal 

will be an eigenstate of 3 , the projection of angular 

.z 

momentum. 

The basis chosen in different shells are given 
below: - 

sd- shells- 0d^ /2 , 1 /2 ; 0d 3 /2 5 i/2’ 1s l/2,1/2’ 0d 5/2,3/2 ? 
0d 3 / 2 , 3 / 2 ’ 0d 5 / 2 , 5/2 

We have taken experimental single -particle energy 
of oxygen-1 7 


Orbit 

Energy (Mev) 

0d 5/2 

0 

1s 1/2 

0.87 

0d 3/2 

5*o8 
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f-p shell 

The basis is 


f 7/2, 1/2’ f 5/2,1/2’ p 3/2,1/2’ 

p 1/2,1/2’ f 7/2, 3/2’ %/2,3/2’ 

p 3/2,3/2’ f 7/2,5/2 ; f 5/2,5/2’ 

^7/2 ,7/2 


Single particle energies 

Orbits 

Energy (Mev) 


0f 5/2 

6.25 

- 

1p l/2 

3*9 


1p 3/2 

2.1 


° f 7/2 

0.0 

f-p-g shell: - The b asis is 



0g 9/2,9/2’ 0g 9/2,5/2’ ° f 5/2,5/2’ ° g 9/2,1/2’ ° f 5/2,1/2’ 

1p 3/2,1/2’ ° g 9/2,-3/2 ; °%/2 f - 

-3/2’ 1p 9/2-3/2’ 

0g 9/2-7/2 

Single particle energies:- 

Orbits 

Energies (Mev) 


1 p 3 /2 

0.0 


C\J 

or 

0 

0.78 


1 p 1 / 2 

1 .05 


0g 9/2 

3.25 


III .3 Storage of Clebsch Gordan coefficients:- 

For storing all sorts of Clebsch-Gordan Coefficient, 
we use the indexing procedure , 
n 1 = C^+1/2) 2 + m, +1/2 

n 2 = ^ 2 +1/2)2 + m 2 + 1/2 


( 111 . 2 * 1 ) 



Let us see how n^ goes, 


when, 3 = 1/2 = - 1/2 ^ = 1 

0 1 = 1/2 m 1 = 1/2 Oj = 2 

j-l = 3/2 = - 3/2 Oj = 3 and so on. 

Similarly for 3 . 

The Clebsch-Gordan coefficients are again stored as a 
three-dimensional array of numbers; viz., CLBGD (n^ ,n2 ? J) 
where J is the total angular momentum. The Clebsch-Gordan 
coefficients are generated by a subroutine called THREEJ 


III . 4 - Generation of trial wavefunction to start Hartree-Fock 
iteration. 

We generate the trial wave functions of the EF 
iteration by diagonelising the matrix of a quadrupole 
operator, since the deformation of nucleon orbits are 
quadrupole type. In (nljm) harmonic oscillator basis, 
the matrix element of quadrupole operator 



1 6 
5 


2 2 

r y o 


) 


lay b e written a s 

p 1 „ " s \ 


nljm 


q; 


i 16 jr 1 

■ 1 

(23+1) (23 1 +1)5 

/ f i 

| h _J 


t 1 

n 1 




I nl> can be easily evaluated utilising the 
. / 
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property of harmonic oscillator hamiltonian and so we can 
write the matrix element above as the product of t wo 
Clebsch-Gordan coefficients and a constant • The matrix 
to be diagonalised is completed by adding the single- 
particle energies to the diagonal elements. The diagonali- 
sation and finding out the corresponding eigenvectors are 
done by a subroutine called EIGEIT. 

III. 5 The Eartree-Fock iteration. 

The KF iteration is done by comparing the HF energy 
Egp at every stage of iteration. Egj, is calculated with 
the IIF wavefunctions at every stage of iteration. The 
Hartree-Fock energy and the quadrupole deformation parameter 
are computed by a subroutine called QERG. At each stage 

of iteration the summation } C ^ L : , is computed by a 

vr -M 

function RHO. 


III. 6 Degeneracies of HF wavefunction. 

If for a particular occupied orbit ^/y > > 
the orbit |P\^> '• £ y 1 ^/. is also occupied then 

the HF hamiltonian is said to have time reversal symmetry. 
An orbit is degenerate with its time-reversed orbit. 

In even-even nuclei, time-reversal symmetry is 
obtained by requiring that, for each occupied orbit, 

^ 1 Vn £ T. y\ _ _ 






'V'lir vX 


•vi v v -m 


' (III. 6.1) 
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the time -reversed orbit 




TIV 


> c.* . 

i ^ '") Tf' 

"h v- }'m 


(V ) * ^ 


/ ; 

"V, 1 / 


,vy (II I. 6. 2) 


V/. 


is also occupied and orthogonal to 

If a nucleus contains equal number of protons 

and neutrons outside a core then we can fill the same orbit 

J/yWith a neutron and a proton. Then the BF hamilonian 

is said to have isopin invariance. In this case the neutron 

and proton orbits are degenerate. In an even-even nucleus, 

if for each occupied orbit the time reversed orbit is also 

filled then tho HF wavefunctions have' four -fold degeneracy. 

( 6) 

III. 7 Angular Momentum projection;- 

, We'll briefly outline the angular momentum projection 

technique here which is actually much more involved. There 

is well-established technique of angular momentum projection based 

(4,5) 

on the framework of Hartree-Fock Bogoliubov (BFB) theory 
where the orbits are a s signed fractional occupation -and 
probabilities (between 0 and 1) unlike in the case of pure 
Hartree-Fock theory where the occupied orbits have occupation 
probability one and others zero 

The intrinsic deformed state in the HFB theory is 
given b y 

t j -r 



TT (ur f b a ^ b; v, ) IV 


L'm 


where the amplitudes u. 1 * and 5T characterise the 

JL v 

emptiness and fullnes^s respectively of the orbit im. 
"U^-Oand "^^=^we get pure HF case. | 0y> is the true 


When 


vacuum of the nucleus. 
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Iiie creation operators are expanded, in terms of 


those in the shell model spherical well i.e. 

! . V i'' -t- i t ■r i - n i ^ 

~ l_L) ^ 


yh 


A 


n - n 


Mil . 7. i) 

(JGT 7 - 5 > 


The wa-efunction (111.7*1) can he written as 


!Y>, u - ■{ i j, ^cc^)| 

\ 


O'- V- (;■' *( ‘ 


'K K 

uQ,' 


- ' H V* 

cx,|_a A— 


;-*v 


A 


' / 


A-'i '• 


h "' i- 




>-A- ■ i -Ar ) 


K - '' ’ f "U. ^ 

We coulri write this because when the exponential is expanded 

in a power series the extra terms vanish identically. We 

have denoted, here, the quantum numbers ( ■) ) by and 

N is a normalisation constant. 

The shell model hamiltonian is 

H ' I (- OK? r "-< + U <-<PN|Yf.> a^r-Zr^a, 

<X r c^|37^ P v x 

So, we can write the energy of a state with angular 
momentum J in the ground state rotational band as ( c.f. 
eq. (II. 2. 10) 

\ h ?; \ f_o> . 


7 

n, > • ' • - 1 - -^-ae 


r" z * (». (r) 


■riVi 


(jr. 7 C) 


T 






O 


The overlap integrals h(9) and n(0 ) are given by 
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and 



respectively . 



( 111 . 7 . 8 ) 



lwJ ’ M - F-f 1 - - ' (Ii.| o) 

The computation goes as follows. First by carrying out 

1IFB calculation we get the intrinsic state \ 7 ^,^ from 

this we set up the f -matrix. Then we compute the capital F 

matrix. This is done by the subroutine CAPF. The M-matrix 

is computed by the subroutine FSMF. 1 / 1 +M and M/ 1 -Hi (i.e* 

1 ~(1 +M)~*' ) are computed by the subroutine MATIOT. Then 
the overlap integrals and are computed by the 

subroutine, TBM. Finally we get Bj from the expression 


(III. 7. 6 ) 



CHAPTER IV 


RESULTS AND DISCUSSIONS 

IV. 1 Is-Od Shell:- 

The Is-Od has, since long, been considered a region 

, ( 1 ) 

of well-df formed nuclei. As pointed out by Ripka , one 
can assod ate w ell defined intrinsic states with the 
1 s-Od shell nuclei. 

The 1 s-Od shell nuclei have been extensively 
studied in the frame work of the exact shell model in 
the recent past^. It is thus not necessary to examine 
these nuclei in the framework of the phenomenological a 
models such as IBM. 

But, here we have found that the collective 

nature of the low-lying levels in the 1 s-Od shell nuclei 

do indeed depend selectively on the V(J=0+2) part of 

the interaction. From Figs. (1-5) have presented the 

22 2 } i 24- 

low-lying collective levels of nuclei ( He, Ne, Mg, 
2 ^Mg, 3°Si) with J=o+2 interactions only, as well as 
full interactions so that the two spectra can be readily 
compared. Since we are mainly interested in the two-body 
part of the interaction, here and elsewhere, we have 
suppressed the role of single particle energies by making 
them degenerate i.o. making all of them zero. 
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IV. 2 Ip-Of Shell: - 

To show that the emergence of the rotational spectra 
in the presence of the V(J-o+2) interaction ( with the single 
particle energies absent) is not an attribute of our calcu- 
lational framework, but an intrinsic proper by of the effective 
interactions ( for the 1 s-Od shell) , we next present some 
results for the 1p-0f shell nuclei. The spectra resulting 
from the V(J=o+2) interaction in the case of the nuclei 
^Ti, ^Fe, and ^Ni exhibit a complete lack of rotational 
spoctra (Figs. 6-8). 

We may thus conclude that the IBM is likely to prove 
extremely poor phenomenological method for studying the 
low-lying levels in the A=4o-70 mass range. 

IV. 3 of-1p-Og Shell :- 

In Figs. 9-10 we present the results for the nuclei 
^Kr and ^Se in the 0f-1p-0g shell. 

The results indicate that the qualitative features - 
the approximate J(J+1) behaviour-are unaffected in going 
from the f ull interaction to just its n J=0+2 subset, when « 
the single particle energies are assumed to be degenerate. 

We have also done the calculations for the nucleus 
^°Kr (Fig.11). Here we find that the spectrum obtained 
with the V(J=0+2) interaction is no longer rotational for 
J y If • This feature may be related to the fact that the 
nucleus 8o Kr is the least deformed of the three nuclei. 
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considered, by us in this region, because of t he approaching 
shell closure at (Z=4-0 ,11=50) • 

It should be mentioned here that the spectra in 
Of-lp-Og shell nuclei have actually been analysed in a 
phenomenological manner, with the help of the Interacting 
Boson Model, iiiis work has been reported by Gelberg and 
U . Kaup ^ recent ly . 
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CHAPTER V 
CONCL USION 

In 1 s-Od and Of -1 p-Og shell nuclei, the structure 
of the low-lying collective levels depends selectively 
only on the J=o and J=2 nucleon pairs. In 0f-1p shell 
nuclei, the low-lying collective spectra are not 
reproduced when one employs the effective interaction 
operating only in the J=0 and J=2 states. 

Since the selective dependence of the collective 
features on only the J=0,2 subset of interaction is an 
important ingradient of the Interacting Boson Model, 
one can infer , that the model is likely to be more 
successful in the 1 s-Qd and 0f-1p~0g shells than in the 
Of -Ip sholl * 
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